Separation of Transitions with Two Quantum Jumps from Cascades 
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We consider the general scenario of an excited level \i) of a quantum system that can decay 
via two channels: (i) via a single-quantum jump to an intermediate, resonant level |m), followed 
by a second single-quantum jump to a final level |/), and (ii) via a two-quantum transition to a 
final level |/). Cascade processes |i) — > |m) — >■ |/) and two-quantum transitions |j) — > |m) — >■ |/) 
compete (in the latter case, \m) can be both a nonresonant as well as a resonant level). General 
expressions are derived within second-order time-dependent perturbation theory, and the cascade 
contribution is identified. When the one-quantum decay rates of the virtual states are included 
into the complex resonance energies that enter the propagator denominator, it is found that the 
second-order decay rate contains the one-quantum decay rate of the initial state as a lower-order 
term. For atomic transitions, this implies that the differential-in-energy two-photon transition rate 
with complex resonance energies in the propagator denominators can be used to good accuracy even 
in the vicinity of resonance poles. 

PACS numbers: 12.20.-m,12.20.Ds,31.30.jc 
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I. INTRODUCTION 

In this article, we consider quite a general problem 
which is illustrated on the basis of the radiative decay of 
excited atomic levels. Let us suppose that an initial state 
\i) of a quantum system can decay into a final state |/) 
via an intermediate, virtual state |m), under the influ- 
ence of an interaction potential V , with relevant matrix 
elements Vfm and Vmi- If all available levels \m) are 
nonresonant, then the decay rate can be computed using 
time-ordered second-order perturbation theory One 
famous example is the decay of the 25 state of hydrogen, 
whose decay to the ground state is dipole-forbidden. Nev- 
ertheless, the main contribution to the 25 state decays 
to the ground state is caused by the very electric dipole 
coupling of the bound electron to the quantized electro- 
magnetic field: the transition proceeds via virtual nP 
levels (n > 2), which are all nonresonant in the nonrela- 
tivistic approximation. One has to formulate the problem 
in second-order as opposed to first-order time-dependent 
perturbation theory. 

Cascade decay accompanying the process \i) \m) — !> 
I/) can proceed when some of the available virtual levels 
\rn) G arc resonant. In that case, the atom may 

first undergo a transition \i) \rn), then \rn) — > |/) (cas- 
cade decay). An example is the decay 35 — ^ riP —5-15 
in atomic hydrogen, where the atom may first radiate a 
photon at the resonant frequency of the 35 — )- 2P tran- 
sition, and then radiate a second photon at the resonant 
frequency of the 2P — > 15 transition. However, the tran- 
sition 35 — >■ nP — > 15 may also proceed via a nonreso- 
nant nP level, in which case it is a true two-photon (two- 
quantum) transition. Indeed, the second-order transition 
amplitude for the sum of the processes 35 — > nP — >■ 15 
(n being summed over) contains both the transition am- 
plitude due to nonresonant virtual states as well as the 
transition amplitude due to resonant intermediate states. 
Strictly speaking, the situation is even a little more com- 
plicated: the electric-dipole coupling of the atom with 



the radiation field couples a state with the atom in the 
35 state to a combined atom+field state with the atom 
in the \nP) state and one photon in the radiation field. 
We can denote this state as |m) = \nP, Ij:-^) for the par- 
ticular transition mentioned. Here, k is the photon wave 
vector, and A is its polarization. Unless simultaneously 
n = 2 and the photon fulfills the resonance condition 
^'35 ~ = Tick, where k = \k\ is the wave num- 
ber of the photon, the intermediate state |m) is non- 
resonant. An intermediate level |m) = |2P, 1^:,^) with 
E^s ~ E2P ^ Tic k' constitutes a nonresonant level even if 
the atomic part of the intermediate state — the 2P level — 
can become resonant. To give another example, an inter- 
mediate state |m) = |4P, 1^_^) with arbitrary k is always 
nonresonant because there is no photon frequency avail- 
able which could turn this level into a resonant state. The 
question then is how to separate the decay through reso- 
nant intermediate states from the decay via nonresonant 
intermediate states. Certainly, it is impossible to do this 
by excluding the |2P) level from the sum over the inter- 
mediate atomic levels, because this level can be both res- 
onant (if the photon frequency in the intermediate state 
is resonant with respect to the 35 — > 2P transition) or 
nonresonant (if the photon frequency in the intermediate 
state is nonresonant with respect to the 35 — ?> 2P transi- 
tion) . The exclusion of the 2P state had been proposed 
in Ref. [2| but has since been scrutinized 0, 1^ . 

Related questions are investigated here in more general 
terms: How can we formulate the problem, within time- 
dependent second-order perturbation theory, so that the 
resonant intermediate levels in the process \i) — > |m) — > 
I/) are separated from the nonresonant levels, and so that 
the cascade contribution due to resonant intermediate 
levels \rn) is clearly identified within the time-dependent 
formalism? In order to answer this question, we first re- 
call that under rather general assumptions about the pro- 
cess, the intermediate states |m) represent a continuum 
of states. This is the case even in transitions of discrete 
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atomic levels because the intermediate states |m) in this 
case are product states of the atom in a discrete state 
and one or more excited modes of the electromagnetic 
field. While the bound states of the atom are discrete, 
the photon modes represent a continuum of energies. In 
particular, the photon wave vector k represents a continu- 
ous variable. A resonant process involves a transition to a 
lower atomic level with a simultaneous emission of a pho- 
ton of the resonant frequency; in that case, the resonant 
state |m) is an eigcnstate of the unperturbed Hamiltonian 
of atom-|-radiation field with exactly the same energy as 
the initial state (the sum of the energies of the lower 
atomic state and of the energy of the radiated photon is 
equal to the energy of the initial atomic state). 

When the decay \i) \m) — >■ |/) can proceed via 
a resonant state |m) which can be reached from \i) via 
a single quantum jump, we have to take into account 
both possibilities: (a) the sequential transition (cascade) 
and (b) the two-quantum transition via the nonresonant 
levels. One possibility to identify the cascade within 
time-dependent perturbation theory is given by the func- 
tional form of its time dependence: for a cascade decay 
\i) — I to) I/), the probability of finding the system 
in the final state |/) is proportional to the square of the 
elapsed time t: The system first has to make a quantum 
jump \i) — > I to), leading to a linear increase (with time) 
of the population of the resonant intermediate level |to). 
The second quantum jump |m) — |/) then leads to a 
quadratic increase of the probability of finding the system 
in state |/) with time. By contrast, the true nonsequen- 
tial two-quantum transition \i) — >■ |m) — > |/) via nonres- 
onant intermediate states leads to a linear increase (with 
time) of the probability of finding the system in state 
I/) with time. Here, we identify, in a general formalism, 
those contributions of the two-quantum transition which 
contribute to the linear behavior (in time), and separate 
them from the (quadratic in time) cascade effect. 

We follow Ref. [l[ in our conventions and proceed as 
follows: First, the basics of a single-quantum transition 
are recalled (Sec.|ll|. We then proceed to the discussion 
of a two-quantum transition without cascades (Sec. 
before including the cascades/resonant levels in Sec. IIVI 
Conclusions are reserved for Sec. |Vl The interaction is 
switched off adiabatically in the distant past and in the 
distant future, but the rate is calculated near t ^ 0. We 
work in natural units (?i = c = eo = l)- 



II. SINGLE-QUANTUM TRANSITION 
A. General formulation 

Following Chap. 5 of Ref. Q, we first consider a single- 
quantum transition \i) to |/), with Cf(t) being the time- 
dependent expansion coefficient of the final-state Hilbert 
vector with respect to the state |/). The interaction is 
adiabatically damped on in the infinite past t — )■ ~oo and 
suppressed by an exponential factor exp{rit), with 77 > 



being an infinitesimal parameter. We then start the time 
evolution with c/(0) = and Ci{0) = 1 (initially, the 
system is in the state \i)). For the complex probability 
amplitude Cf(t) of finding the system in state |/) at time 
t, one finds [see Eq. (5.8.2) of Ref. 



c/(0 = -i J V7,e''*'e'"^'*' dt' 
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Vf^, (1) 



where Vf t is the matrix element of the interaction Hamil- 
tonian V in the Schrodinger picture, i.e., Vfi = {f\V\i). 
Note that there is a somewhat subtle difference between 
the interaction Hamiltonian V in the Schrodinger picture, 
and the interaction Hamiltonian cxp{iHQt)V exp{—iHot) 
in the interaction picture, because in the latter case, ma- 
trix elements of V acquire a time dependence. This time 
dependence is explicitly written out in the term Vfi e"^^'* 
in Eq. (P). 

In the case of an electric-dipole transition in an atom, 
V is the coupling of the bound electron to the quantized 
radiation field. The expression ujfi = Ef ~ Ei is the 
energy difference of the initial and final state of the tran- 
sition with respect to the unperturbed Hamiltonian 
of the system. In the case of an electric-dipole transition 
in an atom, Hq is the sum of the unperturbed Hamilto- 
nian of the atom and of the electromagnetic Hamiltonian 
counting the modes of the radiation field. From Eq. ([T|), 
we find |c/(i)|2 = e^''* |V/i|V('^/» + V^)- Differentiating 
this expression with respect to time, we obtain 



With the identification [see Eq. (5.8.5) of Ref. Q 
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we obtain in the limit 77 —> 0^ 



.(1) 



dt 



|c/(i)P 



77 0^ 



27r\Vf,\^S{L0f,), 



(2) 



(3) 



(4) 



where by definition, F^^'' is the decay rate associated 
with the transition \i) — s- |/) via a single quantum jump 
(we reemphasize that the time derivative is taken at 
t = 0). This result is known as Fermi's golden rule [see 
Eqs. (5.6.35) and (5.8.6) of 

One might wonder why the Dirac S persists in the final 
result, although Fermi's Golden Rule is known to be di- 
rectly applicable to experimentally relevant calculations, 
and an expression containing a Dirac S might otherwise 
be assumed not to be applicable to an experiment. Just 
after Eq. (5.6.35) of Ref. which is equivalent to Eq. (g]) 
in this work, it is stated that the final state must be in- 
tegrated over an (infinitesimal) interval about the final- 
state energy This statement is useful, but it may need a 
more complex explanation for a full elucidation. Indeed, 
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the solution to this question involves two observations: 
(a) that Eq. Q needs to be summed over the state vari- 
ables of the radiated quanta (in the case of an atomic 
transition, photons) in order to make experimentally rel- 
evant predictions, and (b) that the Dirac 6 disappears 
when all possible energies and all possible polarizations 
of the emitted quanta are taken into account in the final 
state. In order to illustrate this aspect, we now discuss 
the application of Eq. Q to an electric dipole transition 
in an atom. 



B. Specialization to an atomic transition 

In the case of an electric-dipole transition of an atom, 
the final state |/) is a product state of the atom in state 
I/a) and one radiated photon jl^:^) in the radiation field. 
In the following, we will write a general product state |/) 
of the system composed of the atom-|-radiation field as 



I/) = IfAj) 



(5) 



where is atomic part of the product state, and |/) is 
the photon part of the product state. The unperturbed 
Hamiltonian of the system is 



Ho^J2 ^f- 1/^) (f^ I + E ^ 4a «£a ^ 



(6) 



fcA 



where the a^^ and at^ are photon annihilation and 
creation operators (here, we work in a representation 
with a finite normalization volume V, i.e., [aj:_^,ai^J = 
Sxx' ) . Eigenstates of the Hamiltonian ([5]) are product 
states of the atom in eigenstate and a Fock state of 
the electromagnetic field such as l^j:^)- 

In the Schrodinger picture, the dipole interaction of an 
electron at point x with the quantized electromagnetic 
field is given by 



V= -ex-E, 



(7) 



kX 

with the unit polarization vectors e^^ and the electric 
field operator E. In atomic physics, one distinguishes 
between the {p ■ A) and {x ■ E) forms of the interaction 
with the electromagnetic field. The former is called the 
velocity gauge because of the appearance of the electron 
momentum in the interaction Hamiltonian. The latter is 
commonly referred to as the length gauge, because the 
electron coordinate x in the interaction Hamiltonian has 
physical dimension of length. In some situations, the 
length gauge is preferable because the interaction is for- 
mulated in terms of physically observable electric field 
strength E instead of the gauge-dependent vector poten- 
tial A (see Refs. [l-Q). AH results presented here are 
given in the length gauge. 



While the initial state of the atom and the final 
state of the atom |/^) are well-defined for an atomic de- 
cay rate, we have to sum over the degrees of freedom of 
the radiated photons in order to obtain the decay rate 
for the one-photon transition — > \fA)- The Dirac 
6 function in Eq. ^ ensures the fulfillment of the res- 
onance condition. For a one-photon final state, we can 
replace 



(8) 



fcA 



for the sum over the photon degrees of freedom of the 
final state. Indeed, the atomic one-photon (I7) decay 
rate for the transition |/a), which we denote as 



r 



(I7) 

/a iA 



is obtained as the sum 



/a iA 



\iAfl) 



(9) 



kX 



kX 



27rY,SiEf,-E,,~k)\{fA,lf:^\V\iA,0)\^ , 



fcA 



where we recall that the sum over k and A transforms 
into an integral in the continuum limit. 



e-^e/ 

fcA ^ 



(27r)3 



(10) 



This integral cancels the Dirac 6. We reemphasize that 
the Dirac S function is eliminated after a summation over 
specific degrees of freedom of the final state, namely, the 
degrees of freedom of the electromagnetic field. 

One might wonder why the single-quantum transitions 
apparently conserve energy according to the above for- 
malism [persistence of the S{ujfi) in Eq. Q], while spon- 
taneous decay of an atomic state always tends to lower 
the energy of the bound electron. The answer is that the 
final state of the process, which is a bound electron in a 
lower state plus a single resonant photon, has the same 
energy as the initial state (electron in the excited state 
and no photon in the radiation field). This is manifest in 
the expression S{ujfi) = 5{Ef^ ~ Ei^ — k) in Eq. 

In view of Eq. ([7]), the transition matrix element 
(/yi, 1^:^ \y\ lA, 0) in Eq. ^ can be written as 



^kX 



-e X 



E 



M,0) = ~e\I^Jj:^ ■ {fA\x\iA) , 

(11) 

where we denote the atomic component of the bra and ket 
vectors by a subscript A. The sum over the photon modes 
in Eq. ([7]) collapses because there is exactly one definite 
photon mode occupied in the state l^;^)■ Summing 
over the available photon modes in the exit channel, we 
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obtain {k 



and obtain 



kX 



fcA 



E 



■47r^afc le^:^ • {/a \x\ m)| 



(27r)3 ^" l^'^'^ 

2a{Ef^ ~ E,^)H^'^^ {/a \x^iA) {lA \x''\ /a) 

(12) 



47r 



^aiEf~E,f \{fA\x\iA)f 



which is the famihar result for a one-photon electric- 
dipole decay rate. The transverse deha function is 
5'^'''^ = S^^ — k'-y /k'^. We denote the Cartesian com- 
ponents of a vector by superscripts. Note, in particular, 
that the sum over the photon modes in Eq. ([7]) is not 
enough in order to calculate the familiar expression for 
the one-photon decay rate; an additional summation over 
final states is necessary. 



III. TWO-QUANTUM TRANSITION 
WITHOUT CASCADES 



A. General Formulation 

In second-order time-dependent perturbation theory, 
the amplitude Cf{t) to find the system in state |/) at 
time t due to the transition \i) — > |m) — > |/) is given by 



cf{t) = dt'e"*'+'"^"*V/„ / dt"e"*"+'"''"*'V, 



which leads to [see Eq. (5.6.37) of Ref. 



{Alf + u%) 



E 



Vfm Vrr, 



(13) 



(14) 



This is a generalization of Eq. ([T]) to second order. When 
no cascades are allowed, we can differentiate with respect 
to time and assume that w,„i is always nonvanish- 
ing. In oder to fix ideas by comparison to a concrete 
example, we recall that in the case of the 25* — ^ 15 two- 
photon transition in atomic hydrogen, the intermediate 
state |m) = 1^^^) has a higher energy than the initial 
state \i) = \2S, 0) where the atom is in the 25 state and 
the electromagnetic field is in the vacuum state |0). No 
cascades are relevant in this case, and Eq. is imme- 
diately applicable. 

We can thus differentiate Eq. ([T^ with respect to time 



.(2) 



d^l^^l 



477 



t=0 



E 



2'!TS{uJfi) 



2 



E 



(15) 



In analogy to the single-quantum transition described by 
Eq. the Dirac S disappears when the final states are 
summed over the experimentally relevant degrees of the 
radiated quanta. We now verify that Eq. ([15]) exactly 
reproduces the known expressions 0, Q for two-photon 
decay rates in atoms. 



B. Specialization to an atomic transition 

For a two-photon transition in an atom, we can write 
the initial state as \i) — |iyi,0), where [ia) is the atomic 
final state, and |0) is the vacuum state of the electromag- 
netic field. The intermediate state is \m) = |myi,l^^), 
where the atom is in state Itoa), and the electromag- 
netic field is in the one-photon Fock state |lg;y). The 

final state is |/) = |/^, I^^Ai' IfcaAj)' '^here the \ki\ and 
Xi are the wave vectors and polarizations of the two ra- 
diated photons (i = 1,2). 

In specializing Eq. (jlSp to a two-photon transition in 
atoms, we have to take into account a subtlety, which 
we outline in greater detail because it becomes relevant 
for all discussions in the following. Namely, the atomic 



(2) 

decay rate is obtained after summing the rate T j:/ over 
the degrees of freedom of all possible radiated photons. 
Now, if we sum the final states over all fciAi and all 
fc2A2, we count the photons twice, because the Fock state 
1 r \ , 1 r \ ) obtained under the simultaneous exchange 

ki <-> k2 and Ai ^ A2 is identical to the original state 



-'-fciA 



At \ )■ Hence, 

' K2A2' ' 



.(27) 



Ja iA 



(16) 



fci Ai /C2A2 



The factor 1/2 is discussed after Eq. (5.108) on p. 169 of 
the quantum field theory textbook |lO| and in the text 
preceding Eq. (3.316) of the textbook [ll| . 

With reference to Eq. , we now turn our attention 
to the two-quantum decay rate (without cascades). Here, 
two quantum paths are possible which must be added 
coherently. These correspond to a different time ordering 
for the emissions of the photons with photon wave vector 
ki and polarization Xi (i = 1,2). Summing over the final- 
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state photon polarizations, the result then is 



-(27) 

fAiA 



felAi fc2A2 



kiXi /C2A2 

4 h k2 



E 



Vfm Vrr, 



E 



e' S(Ef ~E,-ki- fca) 



^feiAi ■ (/-^ 1^1 rriA)) (%A2 ■ ("^-4 \x\ tA)) 



EniA ^ ^lA + ^2 

(%2A2 • 1^1 "^-4)) (cfc^;,^ • (m^ |a-| i^)^ 



(17) 



Separating angular and radial variables for the photon 
energies, we finally obtain the following known result ^ 
in the continuum limit [see Eq. (jlOp j: 



-(27) 



4 a' 



fA^A-^ J ^\kk\E,,-E.,,-kf 


) {niA \x^\iA) 



E 



{Ja \x^\mA) {niA \x^tj 
EmA ~ + k 



{Ja \x^\mA) {niA \x^\ lA 



y 



~ EfA - k 



(18) 



The integration over k extends over the allowed frequency 
range for a two-photon transition [l^j ■ The subtlety with 
respect to the counting of photon modes illustrates that 
Eq. (|15p cannot be applied to atomic transitions without 
a proper interpretation of all physical quantities involved. 



in the denominator of the expression on the right-hand 
side of Eq. p^ . Again, we have already incorporated 
the infinitesimal imaginary part in Eq. ()19p . Here, we 
extend the discussion beyond that in Ref. [l| and analyze 
the resonant and nonresonant levels separately. 

We now assume that some of the intermediate states 
of the system are close in energy to the initial state of the 
process, i.e. that there exist states |m) with Em = Ei. 
We recall that £Wr here represents the total energy if the 
system. In the case of an atomic transition, this would 
be the sum of the energy of the intermediate atomic level 
and of energy of the photons radiated. In order to ana- 
lyze this process, we restrict, in Eq. (|19p . the sum over 
intermediate states to the resonant states to. Then, 



|c/(t)|^ 



E 



(21) 



where uJmi tends to zero. The cascade contribution asso- 
ciated with the resonant levels |to) needs to be differen- 
tiated twice with respect to the time. We obtain 



dt^"^^" 



16772 



E 



Vfm Vwi 



Att'^ ^ \ VfmVmi\^ S{ujfm)S{u}fi) 

m 

^ 47r2 ^ \VfmVmi\'^ S{ujfm) S{uJmi) 

rn 

(1) T-(l) - 



v^p(i) p(i) ^ r 

/ y ^ fm^mi — '^f 



(22) 



IV. TWO-QUANTUM TRANSITION WITH 
CASCADES 



A. General formulation 

We return once more to Eq. ([T?|) which gives the result 
for the two-photon decay rate [see also Eq. (5.6.37) of 
Ref. i]. 



\cf{t)\' 



E 



Vfm Vrr, 



IT] 



(19) 



In the text directly following Eq. (5.6.37) of Ref. Q, it 
is stated that the best way to deal with the situation of 
a resonant intermediate state with ujmi ~ is to use an 
adiabatic turn-on of the perturbation that leads to the 
transition. We have already incorporated this adiabatic 
turn-on into Eq. (fT3|) . It is also stated in Eq. (5.6.38) of 
Ref. [l[ that the turn-on amounts to the replacement 



UJ„ 



177 



(20) 



In the last step, we define the expression C/; as the rel- 
evant cascade term which we evaluate for atomic transi- 
tions in Sec. IIVBI below. We also assume that interfer- 
ence terms among the different resonant levels |to) van- 
ish. 

Equation ((22|) is just the expected result: the level |i} 
feeds the resonant intermediate levels |to) with a time 
dependence F^i] t, and the resonant intermediate levels, 
in turn, feed the final state population as 



(23) 

The necessity of the sum over to is also clear, because all 
intermediate resonant levels have to be included. 

Now that we have treated the resonant levels sepa- 
rately, we have to subtract them from the remaining 
expression. We thereby obtain a modified probability 



6 



|c/(t)|^ of finding the system in state |/), 



for every term in the sum. We then obtain 



\cfitW = \cfit)\'-\cf{t)\ 



{Aif + u)^) 



E 



E 



i?7 



(24) 



S" = — 27r S^LUfi) 



Vfm V-fn 



fn ^ 2 /m 



27r(5(w/i) ^ 



mi ' \2fm 



One might think that the subtraction term (second term 
on the right-hand side of the above equation) would im- 
ply, e.g., the subtraction of the intermediate 2P state in 
the two-photon decay of the 35 state of hydrogen. How- 
ever, that is not the case. The intermediate states are 
quantum states of the coupled system of atom-fradiation 
field. As already outlined in Sec. U the product state 
composed of the 2P level and a resonant photon would 
qualify as a resonant state |m), but a 2P state with a 
slightly off-resonant photon would not constitute a res- 
onant intermediate state. Therefore, the 2P state may 
not be taken out of the sum over the atomic-state com- 
ponents of the virtual states. The time derivative of the 
subtracted expression |c/P is 



2T:5{ujfi) 



2 k^%\VfTnV^, 



An^Siujfi) E^nr 1^/™^*!^ ^i^rnt) 



m ^ fm 



m ^ fm 



(T) —S{^fm) 



= - 2TTS{uJrn:i) ^ \Vmif = ^ E 



(1) 

mi ' 



(28) 



d^l^^l 



-0 (V+^?,) 
477 



E 



Vfm Vrr, 



E 



'mi 1^ 
VfmVj 



UJy, 



vq 



(25) 



The second term on the right-hand side of Eq. is 
divergent in the limit 77 —i- 0+ and ujmi — 0. Wc cannot 
proceed without giving a physical interpretation to the 
adiabatic parameter 77. First, in the subtraction term 



S=- 



(2r/)2 



- UJ 



y-' yfm Vr 



m 



oj- 



ml 



IT] 



(26) 



we carry out the limit 7y — > 0+ in the prcfactor; this 
leads to a Dirac S. Then, for the sum over rn, we match 
the adiabatic parameter with the imaginary part of the 
interaction Hamiltonian. The in and out states \i) and 
/) arc assumed to be asymptotic, stable states in the 
infinite past and future within the context of adiabatic 
perturbation theory. Adiabatically, we therefore switch 
on only the virtual intermediate states. We should thus 
replace 



V ^ 2^ fm 



(27) 



In going from the fourth to the fifth line of the above 
equation, we have neglected interference terms. This 
deserves some comments, which we give by way of ex- 
ample. Let us consider a situation with an initial 45* 
state without any photons, and resonant 2P and 3P vir- 
tual states (each endowed with a single resonant photon) , 
and a IS* final state (with two resonant photons). A con- 
ceivable 2P-3P interference term would necessitate the 
final states |/) to be equivalent in regards to both their 
atomic components as well as electromagnetic-field com- 
ponents. However, because the emitted resonant photons 
for 45* — J> 3P — ?• 15* have different energy as compared to 
4S' — > 2P — 15, the interference term vanishes. 

The derivation (|28|) clarifies that the subtraction term 
is nothing but the sum of the one-quantum decay rates 
of the initial state to all accessible resonant intermediate 
states. The result coincides with the lower-order sub- 
traction term found in Ref. [l3 for the two-photon decay 
rate, but the above derivation is much more general. It 
means that under this regularization, the two-quantum 
correction to the decay rate is obtained as 



r 



(2) 
/J 



2TT5{ujfi) 



t=0 



E- 

m. 



Vfm V„. 



2 fm 



Er^ 



(1) 

mi 



(29) 



where we introduce the overlining in order to differentiate 
v'^Jl from V^Jl . The result (f29| is well defined and gauge 
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invariant [l^ . We also note that the total (one-quantum 
plus two-quantum) decay rate of level \i) thus is 



photons, we obtain 



2t: S{ujfi) 



E 



^f7n ^mi 



m ^mi 2 fm 



(30) 



This result states, in general terms, that the expression 
for the two-quantum decay rate, in the presence of al- 
lowed cascade transitions and with propagator denomi- 
nators regularized by the total one-quantum decay rate, 
has to be interpreted as a one-|-two quantum decay rate. 



B. Specialization to an atomic transition 



In view of the result given in Eq. (jl8p , it is immediately 
clear how to apply Eq. (j30p to two-photon transitions 
in atoms. Namely, when Eq. ((30|) is evaluated for two- 
photon transitions, the correct result is obtained when 
the virtual-state energies in formula (|18p are regularized 
by their total one-photon decay widths [see Eq. (1) of 
Ref. ^2 for a concrete example and extensive further 
discussion in Refs. [T2l - [l5j ]. For an atomic two-photon 

transition, F^^^"* as written in Eq. (|29l) coincides with 
the imaginary part of the two-loop self-ener gy due to 
cut diagrams with two-photon emission Il6l | , with 
the photons fulfilling the two-photon resonance condition 



It has been shown in Refs. [l^jliBl that F^^' is of the or- 
der of a'^{Za)^ in units of the electron rest mass energy, 
and is thus of the same order as the result obtained for 
two- photon transitions without cascades 0]. It is there- 

(2) 

fore appropriate to refer to F^/ as a two-photon correc- 
tion to the decay rate of an initial state which otherwise 
decays via one-photon decay. For completeness, we note 
that the two terms on the right-hand side of Eq. (|29|) are 
both of order a^(Za)*, but their difference is of order 
a^{Za)^ and thus smaller by two orders of Za. 

The only calculation remaining concerns the verifi- 
cation of the fact that Eq. (|22|) reproduces the prod- 
uct of one-photon decay rates for the cascade process 
m) — > IJa)- We use Eq. (fT6)) in order to sum 

over the two-photon final states and Eq. (fTO]) in order 
to proceed to the continuum limit. Summing the cas- 
cade term over the degrees of freedom of the emitted 



(31) 



= 1 V V Vf^IIf^I^ 

2 / ^ / ^ / ^ fm rm 
fciAifeaAa ™ 

= ^ E _E E E 2'^'^(^/- + + ^2 - - k) 

X 2TT5{Ejn^ - E,^ - k) 

2 



ruA, l,r 



feA 



~ex ■ E 



The summation over k\ is over both polarizations A and 
over an energy interval for k = \k\ which contains the 
resonance frequency of the intermediate atomic resonant 
state After performing the sum over kX and going 

to the continuum limit with the help of Eq. (|10p . we 
obtain 

Cf... = ^'I1J2 _E (2^)' -^(^1 - (^/- - ^-)) 

fciAi fcaAa 



x5(fc2_(£L_^ _^, J) 



fcl 



h |/_ 

X ( TO/1 

2V W 



«-'2 A2 



2V 



fA 



"^fciAi 



TOA ) 



Y.{'^iEf,~E^,f\{fA\x\^A)? 



X Y (^mA - E^J,f \{mA \x\ iA)\^ 



Y^p(l7)_ p(l7) 
/ ^ fAiriA ruAtA ' 



(32) 



This result confirms that the cascade terms can indeed 
be written as the product of atomic one-photon decay 
rates. 



V. CONCLUSIONS 

In this paper, we have reviewed the formulation of 
a transition with a single quantum jump within time- 
dependent perturbation theory (see Sec. Ill Ap . The re- 
sult, which is Fermi's golden rule [see Eq. ([4])], is evalu- 
ated for an atomic dipolc transition in Sec. IIIBI We find 
that the familiar result for the one-photon decay rate 
[see Eq. ((T^ ] is obtained after a summation/integration 
over the degrees of freedom of the emitted photon in 
the continuum limit, as given in Eq. (jlOp . The gen- 
eral formulation of a transition with two quantum jumps 
is carried out in Sec. IIII Al within second-order time- 
dependent perturbation theory. The result for the two- 
quantum decay rate F^'^^ as given in Eq. (jlSp is valid if 
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there are no resonant intermediate states through which 
a cascade decay could possibly proceed. The specializa- 
tion to a transition with two quantum jumps is carried 
out in Sec. lIIIBl where it is shown that a summation over 
the two-photon final states of the process [see Eq. ([T6| ] 
yields the familiar result ([T8| for a two-photon transition 
rate in a hydrogenlike ion (such as the 2S-1S decay). 
Cascade contributions are analyzed in Sec. IIV Al 

By isolating the resonant states within the intermedi- 
ate states of the process, we obtain the cascade contribu- 
tion (1231) after differentiating the probability of finding 
the system in the final state twice with respect to the 
elapsed time. Subtracting the contribution of the reso- 
nant states within the first time derivative of said proba- 
bility, we obtain an expression for the two-quantum cor- 
rection to the decay rate of a system which can simulta- 
neously decay via cascades and two-quantum transitions. 
If we use the identification (|27p for the regularization pa- 
rameter r] that parameterizes the width of the interme- 
diate states, then the effect of subtracting the resonant 
states from the two-quantum decay rate is equivalent to 
the subtraction of the total one-photon decay rate of the 
initial state [see Eq. ((29)) ]. The difference of the formal 
two-quantum decay rate (with the propagator denomi- 
nators of the resonant states regularized by their widths) 
and the one-photon decay rate therefore constitutes the 

(2) 

two-quantum correction F^^ to the decay rate, as noted 
in Eq. 

The formal two-quantum decay rate (with the propa- 
gator denominators of the resonant states regularized by 
their widths) therefore constitutes more than its name 
might suggest: namely, according to Eq. ([5(I)) . it is the 
sum of the one-quantum decay rate and of the two- 
quantum correction and therefore constitutes, in some 
sense, a one-|-two quanta decay rate. This finding pro- 
vides a formal and gen eral justification for the obser- 
vation made in Ref. |12| : namely, that the formal two- 
photon decay rate of, say, a 45 state in a liydrogenlike 
ion contains the one-photon decay rate of 45 as a lower- 
order term which needs to be subtracted in order to ob- 
tain the pure two-photon correction. It is instructive to 
also remember that competing one- and two-photon de- 
cays are not restricted to atomic hydrogen, but also occur 
in other atomic systems of fundamental importance like 
rubidium [T7| . 

Finally, in Sec. IIVBI we analyze the cascade contri- 
bution for atomic transitions and show that under the 
appropriate normalization of the two-photon final state 
[see Eqs. (fTB)) and (PT|) ]. the familiar result is obtained 
[see Eq. (|3H)]. The results reported here have an interest- 
ing consequence, because they imply that the one+two- 
photon decay rate (PH]) gives the correct differential-in- 
energy rate coefficients for recombination codes [H, , 
where the spectrum of emitted photons from both two- 
photon and cascade transitions is needed over the entire 
range of resonant and nonresonant frequencies. In partic- 
ular, the corresponding expression (j30|) can be used even 
in the vicinity of resonant bound-state poles, provided 



these are regularized by their respective decay widths. 

One particular limitation of the treatment discussed 
here is immediately obvious. Adiabatic perturbation the- 
ory does not make any statement about the time points 
t' and t" where the two photons are emitted. Thus, the 
total regularized two-photon decay rate as described by 
Eq. (PU)) contains both correlated emission (with a small 
time difference \t'—t"\) and also sequential emission (with 
a large time difference \t' — t"\). Since it is only the pho- 
ton spectrum, not the correlation [T8| . that matters for 
cosmological recombination, we can use the regularized 
two-photon decay rate (j30p . differential in energy, with 
good accuracy for recombination codes [l^, [l^ which 
need the differential-in-energy spectrum of the emitted 
photons as input. In particular, this means that it is 
not necessary to distinguish specific contributions to the 
complete rate (jBO]) : this rate contains both cascade pho- 
tons and correlated two-photon processes in a natural 
way. Therefore, it is understandable that the authors 
of Ref. [l^ could not give a unique value to their pa- 
rameter Auj; this parameter was introduced in order to 
distinguish between cascade photons and two-photon de- 
cays. If one would like to make a more refined distinction 
between cascade photons and correlated two-photon pro- 
cesses, then one has to go beyond adiabatic perturbation 
theory and analyze the dependence of the process on the 
emission times t' and t" , including loss of correlation as 
a function of |i' — t"\, which may be process- as well as 
environment-dependent (e.g., there may be a dependence 
on the average mean free path of the atoms in their en- 
vironment). 

Our article illustrates both the usefulness but also the 
limitations of adiabatic perturbation theory. Namely, if 
we apply the regularization (P7| consistently, to both 
the first term as well as the second term on the right- 
hand side of ((25)) . then the subtraction term 5 defined 
in Eq. (j26p attains a finite value and can be evaluated 
in closed form [see Eq. ([^5]) ]. Furthermore, as shown 
in Ref. [12|, if the regularization (P7| is applied to a 



two-photon transition in atoms, then there are signifi- 
cant cancellations between the two terms on the right- 
hand side of ((^ . which are both of order a^{Za)'^, but 
their difference is of order a^(Za)^, where Z is the nu- 
clear charge number, and a is the fine-structure constant. 
Two-photon decay rates are of order a^{Za)^. As evident 
from Eq. (11) of Ref. [ll] and from Eq. (^5)) in the current 
work, the cancellation of the lower-order terms depends 
on the particular choice of the regularization. Within adi- 
abatic perturbation theory, the regularization (j27p thus 
appears to be the only one which leads to a consistent 
removal of the infinities that plague the two-quantum de- 
cay rate in the presence of allowed cascade transitions. 
Therefore, our article offers — for the first time in the liter- 
ature, to the best of our knowledge — a connection of the 
adiabatic parameter ry used in time-dependent perturba- 
tion theory to a physical concept, namely, the lifetime of 
virtual intermediate states. 

An interesting connection to the theory of energy shifts 
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of atomic levels can be drawn. Low |20l - |22| observed that 
the calculation of energy shifts of excited states of hy- 
drogenlikc ions becomes problematic at order o?(Zaf' ^ 
due to interference effects of the resonance line shapes of 
atomic levels of different principal quantum number. It 
has been argued that at order o^iZaf" ^ two- loop energy 
shifts of excited states cannot be uniquely associated any 
more with a particular atomic level, due to the predic- 
tive limits of adiabatic perturbation theory (20l - |22| . The 
decay rate at order o?'{Zaf' constitutes the imaginary 
part of the energy shift of that same order. It is thus not 
surprising that its calculation requires considerable effort 
within the formalism of adiabatic perturbation theory. 
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